Symmetron dark energy in laboratory experiments 
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The symmetron scalar field is a matter-coupled dark energy candidate which efi'ectively decouples 
from matter in high-density regions through a symmetry restoration. We consider a previously 
unexplored regime, in which the vacuum mass fi ~ 10~^ eV of the symmetron is near the dark energy 
scale, and the matter coupling parameter A4" ~ 1 TeV is just beyond Standard Model energies. Such 
a field will give rise to a fifth force at submillimeter distances which can be probed by short-range 
gravity experiments. We show that a torsion pendulum experiment such as Eot-Wash can exclude 
symmetrons in this regime for all self-couplings A < 1. 



I. INTRODUCTION 

Evidence for a large-scale acceleration of the cosmic ex- 
pansion is now solid [H-Qj but we have yet to determine 
its cause. Current data are consistent with a "cosmo- 
logical constant" vacuum energy density as well alterna- 
tive explanations known as "dark energy." Qualitatively, 
dark energy may differ from a cosmological constant in 
one of two ways. Its energy density may evolve by a 
factor of order unity at recent times, or it may couple 
to known particles more strongly than gravity, provided 
that the resulting fifth forces are screened locally. 

For scalar fields, the simplest dynamical dark energy 
models, several screening mechanisms have been found 
to evade laboratory and solar system searches for fifth 
forces. In a chameleon model 0-01: a nonlinear self- 
interaction makes the effective mass of the field large 
in high-density environments, reducing the length over 
which the fifth force operates. A Galileon field [8l-[ll| 
has a non-canonical kinetic energy, causing it to decou- 
ple from matter at high densities. The symmetron is 
a canonical scalar whose effective potential is symmetric 
under ^ — > — ^ ^M- regions of low density, the symme- 
try is spontaneously broken, and the effective couplings 
are proportional to the vacuum expectation value (VEV) . 
At high density the symmetry is restored, the VEV be- 
comes zero, and couplings vanish. 

If the density of a laboratory vacuum is low enough for 
the symmetron field to enter its broken-symmetry phase, 
then the field will mediate a fifth force between massive 
objects in that vacuum, which may be probed experi- 
mentally [l3l - [l7| . For a symmetron mass n ^ 10"^ eV 
and a matter coupling energy M ~ 1 TeV, this sym- 
metry breaking will occur at densities p < /i^M^ ~ 
0.1 g/cm"^ and distances p,~^ ~ 0.1 mm readily ac- 
cessible to short-range gravity experiments such as the 
Eot-Wash torsion pendulum |18| . In this work we solve 
the symmetron equations of motion exactly in a pla- 
nar geometry. Using this result, we generalize the one- 
dimensional plane-parallel (IDpp) approximation of [l^ 
to symmetron models, allowing us to estimate torsion 
pendulum constraints. In particular, we show that all 
self-couplings A < 7.5 are excluded for M = 1 TeV and 
fi at the dark energy scale. A very interesting region of 



parameter space is accessible to existing experiments. 

This article is organized as follows. In Sec. [ll] we de- 
scribe the symmetron mechanism and apply it to one- 
dimensional matter configurations. Section IIIII uses the 
IDpp approximation to estimate constraints on sym- 
metrons, and Sec. IIVI concludes. 



II. SYMMETRON PHENOMENOLOGY 

A. Symmetrons and screening: At low energy 
and in a flat spacetime background, the symmetron ac- 
tion is approximated by 



5* = 



(1) 



where p is the matter density, A is the dimensionless self- 
coupling, and both the coupling energy M and the vac- 
uum mass p have units of energy. The vacuum energy 
resulting from ([T]) by itself cannot lead to the observed 
acceleration, so we must add a constant term M\, where 
A^A = 2.4 X 10-3 eV is the dark energy scale; however, 
such a constant is invisible in laboratory experiments. 

Previous constraints have considered p ^ lO^M^/Afpi 
and M around the GUT scale, leading to unscreened fifth 
forces on cosmological scales [l^, • Here we are inter- 
ested in /i ^ AfA and M ~ 1 TeV, an energy possibly 
associated with physics beyond the Standard Model. We 
will see that such scales result in new effects in laboratory 
experiments. The symmetron equation of motion is 

□</) = Kff,0(</',p), ^cff = - m') ^ + (2) 
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where the subscript ", 0" denotes a partial derivative with 
respect to 0. At low densities p/M"^ ^ p^ , V^ff has a local 
maximum at (/) = as well as two minima (j) = ±(/)br with 



mb, 



V) = P^/2. (3) 



The Z2 symmetry — >■ — is spontaneously broken as the 
field chooses one of these minima. Henceforth we assume 
= 4'hT in the broken-symmetry phase. Meanwhile, at 
high densities, the mass-squared term ttiq = pjM'^ — p^ is 
positive, and (/) = is the only minimum of the potential. 
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Consider an object with p/M"^ ^ p? at rest in a vac- 
uum. If the object is sufficiently smah, we may linearize 
about the VEV, (p = 0br + <5(?!), reducing ^ to V^5(j) w 
p(f>hi/M'^ + 2fi^S(p. For small fi this looks like Poisson's 
equation = p/ (2Afpi) for the gravitational potential 
^, so we have Stp ~ 2MpjAf^^0br5'- Evidently S(j) couples 
linearly to the density with an effective matter coupling 
/3br = (t>h,Mpi/]VP = 6.0A-i/2(^/10-3 eV)(Af/TcV)-2 
which is much smaller than the apparent coupling /3 = 
Mpi/M = 2 X 10i5(M/TeV)-^ The symmetron in the 
linear regime behaves like a Yukawa scalar with mass 
mbr and force strength a = 2/3^^,. Even this attenuated 
coupling /3br is highly constrained by short-range gravity 
experiments, which exclude a > O.I for m^^ > O.I mm. 
For parameter values of interest here, symmetrons which 
remain in the linear regime are excluded except at the 
highest M and A, so we do not consider them further. 

If we consider larger and larger objects, then the 
assumption of linearity is eventually violated; 6(j) be- 
comes large and negative, the "source" VeS,<j> in ^ 
turns off, and the object becomes screened. Since (50 
cannot be less than —(phi, screening is important when 
1*1 ^ M2/(2Mp[). A typical laboratory test mass with 
p = 10 g/cm^, size ~ I cm, and |*| ~ 3 x 10"^'' wih be 
screened for M < 100 TcV. It is this screened, nonlinear 
regime of the symmetron fifth force which we study here. 

At very low M, the symmetron-matter coupling /3 = 
Mpi/M is large, and symmetrons should be visible in 
colliders. Reference [21| computed collider constraints 
on chameleon models and found that the coupling en- 
ergy had to be greater than ~ I TeV. Although a similar 
analysis for symmetrons is beyond the scope of this pa- 
per, it seems unlikely that M < 100 GeV is consistent 
with collider data. Thus we are interested in the range 
100 GeV < M < 100 TeV. 

We have also considered adding a photon coupling 
term (j)'^F^^F^''' /{8M^) to ©, making the symmetron 
accessible to oscillation experiments |22| - [2a |. An anal- 
ysis similar to (27j shows that symmetron-photon os- 
cillation will occur in the broken-symmetry phase 
with an effective coupling /3-^,br = (phrMpi/M^ = 
6A-i/2(^/10"3gV)(M^/TeV)"^ However, oscillation 
experiments probe /^-y^br si 10^" for p < I0~^ eV, which 
requires either a small A strongly excluded by the fifth 
force constraints of Scc. lIIII or a small likely excluded 
by colliders. 

B. One- dimensional planar configurations: 

The symmetron field profile 0(z) can be found exactly 
for a planar gap between two thick planar slabs. Let 
p{z) = for \z\ > Az/2, inside the slabs of matter, 
and p{z) = pv for \z\ < Az/2, the vacuum between the 
slabs. Assume that pv < p^M"^ < p^n so that the matter 
is screened and the vacuum is possibly in the broken- 
symmetry phase. The equation of motion ^ reduces 
to 0,22 = V;ff,0(0, p(z)), where (/),22 ee £(j)/dz^. Using 

Jz^^ ~ 5^(iz)^ integrate the equation of motion 



over any interval (za, zb) with constant p, we have 



t>B,p)-Ve«i(l)A,p), (4) 



with 4>A = (pizA) and (pB = 4>{zb). 

Equation (|4]) is helpful if we choose either za or zb 
such that d(j)/dz — 0. Choosing the interval (0, Az/2) we 
obtain 0s^z/2 = Vof[{(j)s, Pv) — Kff (0g, Pv)- where we have 
defined 4>g = 0(0) to be the field at the center of the gap 
and 0s = 0(A2/2) to be the field on the surface of one of 
the planar slabs. With the interval (Az/2,oo) we have 
'K:ff (0s, Pm), from which we obtain 



AM2 



(5) 



Pm - Pv ° 12(/0m - Pv) ^ 

Next we choose an interval (0,2:) for < z < Az/2 
to determine 0,^(2) /\/2 = -^Vcfi{<i), Pv) - "K^ff (0g, Pv) 
inside the gap. Integrating, we find 0(2) implicitly. 



PvZ 



p. 



rd(f>/V2 



'0 \/Vcs{<j), Pv) - Vcs{(j)g, Pv) 

_ 1 

/ ^2 \ 



2 ' 



(6) 



where pi = p^ - P^IM"^, k\ = 0g/(20g^ " '^g)' and 
F{9, k) = dtj ^\ — k^ sin(t)2 is the eUiptic integral of 
the first kind. We "solve" the gap by guessing 0g, using 
([5]) to find 0s, using ^ with = 0s to find 2, and refin- 
ing our guess 0g until 2 = Az/2. Once the correct value 
of 0g is known. (|5l6p determine everywhere in the gap. 

The quadratic nature of the potential near = im- 
plies that there is a minimum gap size Azmin below which 
0(2) = everywhere. We can sec this by estimating the 
energy change due to pulling the field to some nonzero 
0g inside the gap. The "gradient energy" density associ- 
ated with this change in the field over a distance A2 is of 
order {(pg/Az)'^, while the potential energy density is of 



order 



Thus the net change in energy is negative 



only if A2 > Pv ^- We can find the precise value of A2n 
by considering ^ with = 0s in the limit 0g — > 0, 



Azij^ii^ 

Pv 



sm 

2 



P?M2 



(7) 



Typically pm 3> Pv, pIM^, so the second term inside the 
parenthesis can be neglected, and Azmin ~ t^/ Pv 

Finally, we determine 0(z) inside the slabs of matter. 
Using 111 with interval (2,00) and 2 > A2/2, we find 

A04/12 with ml = 



0,2 = ~\j2Vcs{(p, Pm) = ■ 
Pai/M'^ — p^ . Defining (/^s 



-Vto^02 
= A0s/(I2mg) and integrating. 



0(2) _ 2e-"°(^-'^^/^)(yrTW- 1) 



(8) 



3 



H=3x10"^eV ■ 
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FIG. 1: Estimated torque in the apparatus of [3 for several 
symmetron models. M = 100 GeV for all models shown. 
A = IQ-^ for ^ > 10"^ eV and A = 10"'' for /i < 3 x 10"* eV. 




matter coupling energy M [GeV] 

FIG. 2: Estimated constraints on symmetron dark energy. 
The light-green shaded region is excluded for fj, — 10~^ eV, 
while models in the black shaded region are unscreened. The 
solid red, dashed blue, and dotted purple curves show models 
with n = 3x 10"* eV, 3 x 10"^ eV, and 10"^ eV, respectively 



III. CONSTRAINTS ON SYMMETRONS 

A. One- dimensional plane-parallel approxima- 
tion: The previous section determined the surface field 
(j)s{Az) for arbitrary gap size Az, as well as the field pro- 
file (l){z) inside a planar slab of matter bounding the gap. 
Rather than planar slabs, a torsion pendulum experi- 
ment such as Eot-Wash uses parallel planar disks with 
surface features such as holes. Let the z axis be nor- 
mal to both disks. As a hole on the "source" disk moves 
past another hole on the "test" disk, fifth forces between 
the holes exert torques on the test disk. The IDpp ap- 
proximation [l^ estimates the field at a point {x, y) on 
the surface of each disk by (/)s(Az(x, y)), where Az{x,y) 
is the distance to the nearest point on the opposite disk. 
The field inside the disk is approximated by ^ given the 
surface field. Since (j)s{x,y) will be greater for a region 
directly across from a hole on the opposite disk, there is 
an energy cost to moving the holes on opposite disks out 
of alignment with one another. This change in energy as 
the source rotates is used to predict the torque signal. 

We find the energy per unit area associated with a 
point {x,y) on the surface of a disk by integrating over 
the interior, with the ^-dependence of 4> given by dS]): 

Consider a region on the test disk across from a hole of 
radius rsh on the source disk. The total energy of the 
field in this region is found by integrating ([9]), i?Th = 
/q^'' 2Trr^dr. We can similarly find the energy Egh h^ 
the region of the source disk across from a hole of radius 
rxh on the test disk. Let AEtyi be the energy difference 
between i?Th and the corresponding energy for a region 
on the test disk not overlapping a source disk hole, and 
define A_Bsh analogously for the source disk. Then the 



total energy cost associated with moving the source and 
test holes out of alignment is AEgh + Ai^xh • If the disks 
have -/Vh holes in N-^ rows, then the torque is 

r = ^fscr{AEsh + A^Th) (10) 

where fsa ^ exp(— moztoii) for Zfoii = 10 /im accounts 
for the shielding foil between source and test disks [l9| . 
This IDpp prediction for the signal in a torsion pendulum 
experiment is shown in Figure [1] for several models. 

B. Torsion pendulum constraints: Eot-Wash [l^ 
looked for fifth forces using two molybdenum disks, of 
density ~ 10 g/cm^, in a 10"^ torr vacuum corre- 
sponding to pv ~ 10"^^ g/cm^. Each disk had iVh = 42 
holes in A^i- = 2 rows. Radii of the source and test disk 
holes were rgh = 1.6 mm and ttii = 2.4 mm, respectively. 

Here we estimate constraints from an Eot- Wash-like 
experiment which excludes torques (jlOp greater than 
0.01 fN-m. Eot-Wash probed fifth forces over a range of 
source-test distances 0.05 mm < Azs_t < 10 mm, with 
the strongest constraints typically coming from the short- 
est distances. For fi > 10"^ eV we use fifth force bounds 
at Azg_x = 0.1 mm. At smaller fi, the minimum gap 
size Azjiiin can be several millimeters, so no symmetron 
fifth forces are predicted at the smallest Azs-t- Thus we 
use Azs-T = 2 mm for constraints on /i = 3 x 10~^ cV. 

For ^ = A/a and M = 1 TeV, we find that A < 7.5 
is excluded. Figure [3] estimates constraints on a wide 
range of symmetron models. The IDpp calculation likely 
underestimate actual constraints by a factor of 2-3, since 
only one Azs_t is used for each model, and since (|3 
only counts the energy of the field inside the disks [l9| . 
However, we do not attempt to correct for this here. 

At ^ < 10"^ eV, Azmin is too large to be probed 
by Eot-Wash, though constraints from larger-scale ex- 
periments apply. For p, ^ 10"^ eV, matter becomes 
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unscreened at all M of interest. At low M the Comp- 
ton wavelength is small, and screening by the shielding 
foil [i^ weakens constraints. Meanwhile, at large Af, 
Pm < tJ.^M^ and the gravitational potential of each disk 
*disk < M2/(2Af|i). Determining precise constraints in 
this linear regime is beyond the scope of this work, but 
such models tend to have large fifth forces. 

IV. CONCLUSION 

Symmetron dark energy is an intriguing new model in 
which fifth forces between massive objects are screened 
through a symmetry restoration at high densities. We 
have considered a previously unexplored parameter re- 
gion in which the symmetron self-coupling A ~ 1, the 
vacuum mas fi ^ M\ so that only one small scale is 
necessary, and the matter coupling energy Af ~ 1 TeV 
could be associated with new physics beyond the Stan- 
dard Model. Since symmetrons in this regime give rise 
to fifth forces at distances ~ 0.1 mm, they can be tested 
by current submillimeter fifth force experiments. 

We have found the exact field profile dSMI) in a planar 
geometry. Using the IDpp approximation, we have pre- 
dicted the torque signal due to the symmetron in an Eot- 
Wash-like torsion pendulum experiment. Figure [1] shows 



this prediction for several models and experimental con- 
figurations, and the resulting constraints are shown in 
Fig. [2] for a large range in parameter space. Specifically, 
we have shown that A < 7.5 is excluded for /i = AI\ 
and M = 1 TeV, demonstrating the power of existing ex- 
periments to probe symmetron dark energy. These con- 
straints may be improved in the near future by using a 
numerical computation such as that of [l5| and analyzing 
data from the next-generation Eot-Wash Experiment. 
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